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A four-body calculation of the pnAA bound state, AA H, is performed using the stochastic varia¬ 
tional method and phenomenological potentials. The NN, AN, and AA potentials are taken from 
a recent Letter by Filikhin and Gal, PRL89, 172502 (2002). Although their Faddeev-Yakubovsky 
calculation found no bound-state solution over a wide range of AA interaction strengths, the present 
variational calculation gives a bound-state energy, which is clearly lower than the A H +A threshold, 
even for a weak AA interaction strength deduced from a recent experimental Baa^aH 13 ) value. 
The binding energies obtained are close to, and slightly larger than, the values obtained from the 
three-body dAA model in the Letter. 

PACS numbers: 21.80.+a, 21.45.+V, 21.10.Dr, 13.75.Ev 


In a recent Physical Review Letter [j], Filikhin and Gal 
(FG) described systematic Faddeev-Yakubovsky (FY) 
calculations for the mass number A = 4, strangeness 
S = —2 problem, in which they searched for a particle- 
stable bound-state of AA H. They did not obtain a bound- 
state solution, even for a strongly attractive AA interac¬ 
tion, the scattering length of which is about oaa ~ —3 
fm. On the other hand, they also studied the same system 
by using a three-body dAA model, where the Ad inter¬ 
action was constructed to reproduce the low-energy pa¬ 
rameters of a A pn Faddeev calculation for both the spin- 
doublet and quartet states. In contrast with the four- 
body pnAA calculation that produced no bound state, 
the three-body dAA model produced a particle-stable 
bound-state. One may think that this incompatibility 
raises an interesting problem concerning “the formal re¬ 
lationship between these four-body and three-body models 
which do not share a common hamiltonian .” [l] However, 
we are doubtful that there is really no bound state in the 
four-body pnAA calculation. 

A recent experimental report [§] on the observation of 
AA He in the KEK-E373 hybrid emulsion experiment has 
had a significant impact on strangeness nuclear physics. 
The Nagara event provides unambiguous identification of 
AA He production, and suggests that the AA interaction 
strength is rather weaker than that expected from an 
older experiment Q. 

Before the publication of the Nagara event, we had al¬ 
ready attempted to search for AA H theoretically by per¬ 
forming a complete four-body calculation using a varia¬ 
tional method Si The AA interaction used in those 
studies was strongly attractive with a scattering length of 
oaa ~ —3 fm. We concluded that AA H is particle stable 
provided that the AA interaction is so strong. 

The variational calculation gives an upper bound on 
the energy eigenvalue. In contrast with the FY formula¬ 
tion, although the variational calculation does not nec¬ 
essarily take the exact boundary condition into account, 
the variational principle guarantees that the energy ob¬ 


tained comes close to the exact value from above as the 
trial function is improved. Therefore, starting from an 
identical hamiltonian for the four-body system, if the 
bound-state solution is obtained in a variational calcula¬ 
tion, the exact eigenenergy must be lower than that and 
the FY calculation should achieve this kind of solution. 

In the calculation of this four-body system, determin¬ 
ing the AN interaction is very important. Particularly, 
the strength in the 3 S'i channel of the AN interaction is 
crucial, as well as the strength of the AA interaction, in 
determining whether AA H is particle-stable. 

The purpose of this article is twofold: One is to ex¬ 
amine the recent result of the four-body calculation for 
pnAA by FG. Our four-body calculation gives quite a dif¬ 
ferent result from that of FG, and we discuss the struc¬ 
tural aspects of AA H as a four-body system. Another 
purpose is to clarify the importance of the choice of the 
AN potential in searching for AA H. 

In Ref. Q], FG used phenomenological NN, AN, and 
A A potentials, which have functional forms of a three- 
range Gaussian. The NN potential utilized in the pn 
spin-triplet channel is consistent with the 2 H binding en¬ 
ergy, and the AN potential is parametrized by fitting 
the low-energy scattering parameters for the Nijmegen 
soft-core 97f (or 97e) potential. For the AA interaction, 
since there is no direct information from experiments in 
free space, FG used various parameter sets. A promising 
one, deduced by reproducing the experimental binding 
energy of AA He 0 from an a + 2A three-body model, is 
weakly attractive with a scattering length of oaa = —0.77 
fm. For all of these interactions, the strength and range 
parameters were determined so as to be appropriate for 
5-wave interactions. We thus assume that these interac¬ 
tions are valid only for the even-partial wave component 
of the baryon-baryon interaction in the three- and four- 
body systems. 

For systematic calculations of 2 H, A H, A H*, and AA H, 
we use the same sets of NN, AN, and AA interactions 
as FG used. The Set A AN potential from Ref. 0], 
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FIG. 1: x So and 3 Si phase shifts of AN scattering as a func¬ 
tion of the A momentum, p A . The solid lines are obtained 
from the Set A potential in Ref. [§), the dashed lines from 
NSC97f(FG) in Ref. Q. 


which has a different strength in the 3 Si channel, is also 
used. The parameters of the Set A AN potential were 
determined phenomenologically in order to reproduce the 
A = 3,4 single-A hypernuclei. Figure ^ shows the AN 
S -wave phase shifts. In the low energy region, the 1 S'o 
phase shifts obtained from NSC97f(FG) and from Set A 
are almost identical. On the other hand, the 3 Si interac¬ 
tion of the NSC97f(FG) is more attractive than that of 
Set A. As we show later, both AN potentials reproduce 
the experimental I?a(aH) value, because Ba(aH) is sen¬ 
sitive to the strength of the 1 So AN interaction, while it 
is insensitive to the 3 Si strength of the AN interaction. 
In other words, the experimental information for A = 3 
cannot determine the 3 Si strength of the AN interaction. 
Therefore, the AN interaction used in the calculation of 
AA H has to be tested not only for R A ( A H), but also f° r 
another B\ which is sensitive to the strength of the 3 Si 
AN interaction; for example, one can use _B A ( A H) and 
_Ba( a H*). This is one of the most important points in 
this article, because the calculated Raa value is very sen¬ 
sitive to the choice of the AN interaction, particularly the 
strength in the 3 Si channel. 

In order to check the validity of the choice of the AN 
potential, we calculate A = 3,4 A hypernuclei, using the 
NSC97f(FG) or the Set A AN potential. Only for this 
task, the Minnesota potentially is used for the NN inter¬ 
action. The parameters of the Minnesota potential were 
determined so as to reproduce low-energy NN scattering 
data. The Minnesota potential reproduces reasonably 
well both the binding energies and sizes of few-nucleon 
systems, such as 2 H, 3 H, 3 He, and 1 1 le(j . 

In this work, the few-body calculations of the vari¬ 
ous systems are performed using the stochastic varia¬ 
tional methocMSVM) with correlated Gaussian (CG) ba¬ 
sis functions [g]. The trial function is given by a combi¬ 


nation of basis functions: 

K 

JMTMt = ^2 CkA [^(*1 Ak)XkJMVTM T } ■ (1) 

fc =1 

Here, A is an antisymmetrizer acting on identical 
baryons, x = (aq, ■ • •, xa- i) stands for a set of relative 
coordinates, and XkjM (vtm t ) is the spin (isospin) func¬ 
tion. The spatial part of the trial function, G(x,A ), is 
the CG, which is defined by 

1 A 

G(x, A k ) = exp | — ^Y^akijin - fj) 2 } (2) 

i<j 

1 A_1 

= ex p{-^ ( Ak h Xi ‘ x i}‘ ( 3 ) 

i,j =1 

The (A — 1) x (A — 1) symmetric matrix Ak contains 
A(A — l)/2 independent matrix elements, which char¬ 
acterizes the CG basis and is uniquely determined in 
terms of otkij ■ A set of linear variational parameters 
(ci, • • ■ ,ck ) is determined by using the Ritz variational 
principle. The variational parameters are optimized by 
a stochastic procedure. This is entirely the same as in a 
previous study ||. The reader is referred to Refs. i! 
for details of the calculation. 

Before showing the results of our four-body calcula¬ 
tions for aa H, we report results for the binding energies of 
2 H, a H, and A H* using the same potentials as were used 
by FG. Using the triplet pn(FG) and NSC97f(FG) AN 
potentials, the calculated binding energies were B( 2 H) = 
2.250 MeV, R a (a h ) = 0.237 MeV, and R a (a h *) = 0.010 
MeV. These energies for the three-body systems are con¬ 
sistent with those quoted by FG, although each energy 
is actually slightly larger than that of FG. We think that 
these small discrepancies are due to the s-wave approxi¬ 
mation of the Faddeev calculation. Note that both calcu¬ 
lations for A H* produce a weakly bound state; this means 
that the SVM with CG basis functions and the Faddeev 
calculation with the s-wave approximation do work well 
even for the very weakly bound-state problem. 

According to our previous studies 0. ill 7 AA H should 
have a particle-stable bound state with an isospin of I = 
0, and an angular momentum and parity such that J’ 1 ’ = 
1 + , provided that a strongly attractive AA interaction 
with a scattering length of czaa ~ —3 fm, is used. Using 
such a strong AA interaction, we have obtained a bound- 
state solution for AA H (see Fig. d - The Raa( aa H) values 
(> 1.2 MeV) obtained are more than two-times larger 
than the values obtained in our previous studies (~ 0.5 
MeV) with the Set A AN potential. This is due to the 
difference in the strength of the AN interaction in the 
3 Si channel (see Fig. QJ. 

The four-body calculation using a weaker AA interac¬ 
tion (qaa. = —0.77 fm) is a challenging problem, since the 
three-body dAA model by FG predicts a particle-stable 
bound-state with a very small binding energy (Raa ~ 0.3 
MeV). For such a weakly bound four-body calculation, 
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FIG. 2: Calculated -Baa(aaH) as a function of the scatter¬ 
ing length, asAA- The solid squares were obtained using the 
NSC97f(FG) AN potential and the solid circles by the Set A 
potential. The open squares are the result of the dAA three- 
body model, taken from Ref. jj. The straight lines were 
drawn only for the sake of a guide to the reader. 



Basis Dimension 

FIG. 3: Energy convergence of AA H as a function of the basis 
dimension, K. The interactions are taken from Ref. jj, spin- 
triplet pn, NSC97f(FG) AN, and AA deduced from the recent 
experimental Baa(aaHc). The converged energy is clearly 
lower than the |H + A threshold. 


though the convergence of the energy is rather slow, the 
energy obtained is clearly lower than the A H + A thresh¬ 
old, and we found that the ground state is particle stable 
(see Fig. 0 . This is a genuine four-body calculation, and 
the calculated Saa(aaH) ~ 0.4 MeV is slightly larger 
than that from the dAA three-body calculation by FG, 
as shown in Fig.0 

As can be seen in Fig. 0 the difference in the Baa 
values between the present four-body model and the FG 
three-body model becomes larger as the strength of the 
AA interaction increases. Moreover, the two lines (la¬ 
belled “pnAA” and “dAA”) in Fig. 0 seem to meet each 




(T c ) (Vnn) 

E c 

\/ ( r Aw) a/ ( r Aiv) 

2 h 


18.74 -20.99 

-2.25 

3.85 


aH* 


19.09 -21.20 

-2.12 

3.79 

37.8 

aH 


20.70 -22.30 

-1.59 

3.54 

8.88 

aaH(oaa = 

b- 

N- 

O 

1 

22.28 -23.17 

-0.88 

3.34 

7.92 

aaH(oaa = 

—2.8fm) 

24.73 -24.55 

0.18 

3.08 

4.82 


TABLE I: Energy expectation values of kinetic (T c ) and po¬ 
tential (Vnn) terms, and the sum of these energies ( E c ), for 
the pn subsystem, in units of MeV. The rms distance between 
a proton and a neutron, or between a nucleon and a A, is also 
listed, in units of fm. The spin-triplet pn and NSC97f(FG) 
AN potentials, taken from Ref. []J, were used. 



Ba(X H) 

Ra(!H) 

Ra(1h*) 

NSC97f(FG) 0.24 

2.69 

1.99 

Set A 

0.18 

2.24 

1.14 

Experiment 

0.13 ±0.05 

2.04 ±0.04 

1.00 ±0.04 


TABLE II: A separation energies, given in units of MeV, of 
A = 3, 4 single-A hypernuclei. The Minnesota NN potential 
was used. 


other at the point where oaa = 0 fm. This means that 
the polarization of the pn subsystem is small, and that 
the dAA model is a good approximation if the AA in¬ 
teraction is very weak. The polarization of the deuteron 
subsystem grows as the strength of the AA interaction 
increases. Table [I] lists the energy expectation values for 
the proton and neutron subsystem in each (hyper) nu¬ 
cleus, and also the root-mean-square distances between 
a p and an n, or between a nucleon and a A. Here, T c is 
the kinetic energy of the pn subsystem, which is defined 
by T c = (p 1 —p 2 ) 2 /^niN- The table shows that the in¬ 
fluence of the A particle upon the internal structure of 
the pn subsystem becomes large as the A particle comes 
close to the nucleon. Especially in the case of a strongly 
attractive AA potential, the change in the internal energy 
( E c ) or of the rms distance (\/ (t 2 n1 v)) is significant. 

As can be seen in Fig. 0 the Baa value is sensitive to 
the choice of the AN potential. For the purpose of pre¬ 
dicting whether AA H exists as a particle-stable bound- 
state, the AN potential has to be examined carefully. 
Table 10 compares the Ba values of A = 3,4 hyper¬ 
nuclei. The calculated Ba value of the A = 4 system 
using NSC97f(FG) is larger than that using Set A, or 
larger than the experimental value. Particularly, the 3 S'i 
strength of NSC97f(FG) is apparently too strong to re¬ 
produce the _Ba( a H*) value, though the NSC97f(FG) re¬ 
produces reasonably well the Sa(aH) value. It would, 
therefore, be rash to conclude that A aH has a particle- 
stable bound-state, though the present four-body calcu¬ 
lation with the NSC97f(FG) gives a bound-state solution, 
even for a weaker AA interaction, such as oaa=~ 0.77 fm. 

The present four-body calculation gives quite a differ¬ 
ent result from that of the FY study discussed in Ref. [I[ . 
At present, we have no clear explanation for why the FY 
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search for AA H has not found a bound-state solution. We 
also checked the accuracy of the present variational cal¬ 
culation by examining the virial theorem^; For an exact 
eigenstate of the hamiltonian, H = T + V, we have 


Baa is sensitive to the 3 Si channel of the AN interaction. 
Therefore, a check of the AN potential concerning the 
observed binding energy of only the subsystem, A H (i + ), 
is insufficient. 


(T) = \(W), with W = 

1=1 


( 4 ) 


For the four-body calculation, we obtained the ratio 
2 (T)/(W) = 1.000016 for a weak AA potential (uaa = 
—0.77 fm), and the ratio 0.999962 for a strong AA po¬ 
tential (oaa = —2.8 fm). Therefore, we think that 
the present four-body calculation gives a virtually ex¬ 
act eigenenergy, and that the Baa value obtained by a 
four-body calculation for pnAA should be close to (and 
slightly larger than) the energy given by the dAA three- 
body model. In comparison with the dAA three-body 
model (in Fig. 0 ), the present result seems to be reason¬ 
able, in contrast to that of the FY four-body calculation 
in Ref. |]fl. 

The contribution to the binding energy from the higher 
partial wave components are marginal. The present po¬ 
tentials are all central and have Gaussian form factors. 
This Gaussian radial form (e.g., Ve~ K ^ ri ~ r ^ ) is rewrit¬ 
ten so as to be valid for each angular momentum in terms 
of nonlocal potentials: [§j 


Ve k(T; Tj) __ j drd,r'\S(ri — rj — r'))(S(ri — rj — r)| 
x £*=oy e —£ m l? m (r')lWf). 

( 5 ) 

We also calculated the binding energies in which the NN, 
AN, and AA potentials are restricted to be valid only 
for the l = 0 component. The binding energy calcu¬ 
lated is B( a |H) = 2.382 MeV (a AA = -0.77 fm), or 
B(aaH) = 2.842 MeV (aAA = —2.8 fm). Each energy is 
still lower than the A H+A threshold (For l = 0 truncated 
interactions, we obtained B( A H) = 2.365 MeV). 

We should emphasize that in the study of AA H; the 3 Si 
AN interaction has to be determined very carefully, since 


One might think that the spin-doublet structure of 
A = 4 A hypernuclei is a means of determining the 3 Si 
AN interaction. However, this strategy without any ex¬ 
plicit E admixture would lead us to a serious problem 
concerning the A — 5 anomaly umm- According to 
recent studies, taking account of the explicit E degrees of 
freedom . 1 lA E3- 1 l-lbl. 16(. the AN — EV coupling, plays 
a crucial role in the binding mechanism of s-shell A hyper¬ 
nuclei. In other words, even the spin-doublet structure of 
A = 4 A hypernuclei does not pin down the 3 Si AN inter¬ 
action, and the AN potential used in the study of AA H, 
has to be tested on a complete set of the observed s-shell 
A hypernuclei. Moreover, the a + 2A three-body model 
might be inappropriate for deducing the AA interaction 
in free space from the recent experimental information on 
Baa(aaH 6 )> since AN — "EN coupling plays an important 
role even for ((He, and the rearra ngem ent energy of the 
core nucleus ( 4 He) is rather large [l2,[l2j. Therefore, a 
study aimed at searching for AA H needs not only a four- 
body calculation, but also five-body ( A He) and six-body 
( AA He) calculations. Furthermore, AA — S N coupling 
effects should be explicitly taken into account, because 
the Pauli suppression effect in the 5 channel of a a He is 
appreciably large Q. A theoretical search for AA H is still 
an open subject. 
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